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Abstract. The comparison of quasar absorption spectra with numerically simulated 
spectra from hierarchical cosmological models of structure formation promises to be a 
valuable tool to discriminate among these models. We present simulation results for 
the column density, Doppler b parameter, and optical depth probability distributions 
for five popular cosmological models. 

INTRODUCTION 

A physical picture of the Lya forest in hierarchical cosmologies has recently 
emerged from numerical simulations [1,2] in which the absorbers that give rise to 
low column density lines (Nhi < 10 15 cm" 2 at z ~ 3) are large, unvirialized objects 
with sizes of ~ 100 kpc and densities comparable to the cosmic mean. Since 
the absorbers grow from the primordial density fluctuations through gravitational 
amplification, statistics of the forest may be used to test various models of structure 
formation. We discuss the numerical stability of the statistics against changes in 
simulation box size and spatial resolution in detail elsewhere. [2] We focus here on 
examples from our model comparison study [3] in which statistics of the Lya forest 
are computed in five cosmological models: the standard cold dark matter model 
(SCDM), a flat cold dark matter model with nonvanishing cosmological constant 
(LCDM), a low density cold dark matter model (OCDM), a flat cold dark matter 
model with a tilted power spectrum (TCDM), and a critical model with both 
cold dark matter and two massive neutrinos (CHDM). The initial fluctuations, 
assumed to be Gaussian, are normalized using a 8h -i to agree with the observed 
distribution of clusters of galaxies, although all but SCDM are also consistent with 



the COBE measurements of the cosmic microwave background. By varying a 8h -i 
within a given model (SCDM), we also investigate the dependence of the statistics 
on changes in the fluctuation power spectrum. 

The simulation technique uses a particle-mesh algorithm to follow the dark 
matter and the piecewise parabolic method [4] to simulate gas dynamics. The 
simulation box length is 9.6 Mpc (comoving) with spatial resolution of 18.75/i -1 
(37.5/zT 1 ) kpc for the model (power) comparison studies, respectively. Nonequi- 
librium effects are followed [5] for six particle species (HI, HII, Hel, Hell, Helll, 
and the electron density). We assume a spatially-constant radiation field computed 
from the observed QSO distribution [6] which reionizes the universe around z ~ 6 
and peaks at z ~ 2. Synthetic spectra are generated along 300 random lines of 
sight through the volume, including the effects of peculiar velocity and thermal 
broadening. The spectra are normalized to give a mean optical depth < r >= 0.3 
at z = 3 to agree with observation. We do not include the effects of radiative 
transfer, self-shielding or star formation and so can not address the physics of the 
highest column density absorbers (N HI > 10 16 cm -2 ). 

FIT DEPENDENT STATISTICS 

The synthetic spectra are fit by Voigt profiles at these low column densities to 
obtain column densities and Doppler widths for each line. The slope of the column 
density distribution is insensitive to changes in the size of the simulation volume 
or grid resolution. [2] Fig. 1 confirms analytic work [7] that the slope of the col- 
umn density distribution depends primarily on the power in the model at scales 
~ 100 — 200 kpc and steepens for models with lower power. Each of the five cos- 
mologies in our comparison study agrees with the data at the 3a level, although 
models (SCDM,LCDM,OCDM) with moderate power at these scales are favored. 
The Doppler b parameter distributions are determined not only by thermal broad- 
ening and peculiar velocities of the absorbers, but also by the Hubble expansion 
across their width, and require high spatial resolution to be modeled properly. [2,8] 
The shape of the distribution including the high b tail is well fit by hierarchical 
models. However, as shown in Fig. 2, observations by Kim, et. al [9] yield me- 
dian b parameters significantly higher than predicted for the models favored by the 
column density distributions (LCDM, OCDM, SCDM). 

CONCLUSION 

Although hierarchical cosmologies reproduce the general characteristics of the 
Lya absorption spectra quite well, detailed tests of the models may require new 
methods of analysis for the simulations and observations. Statistics derived directly 
from the flux or optical depth distributions without recourse to any line fitting al- 
gorithm are particularly interesting. [11] For example, the optical depth probability 
distribution, like the column density distribution, is stable to changes in simulation 
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FIGURE 1. HI column density distributions at z = 3 for lines with Nri in the range 
10 12 ' 8 — 10 16 cm -2 . Results for SCDM varying a sh -i (top); results for the five cosmological models 
(SCDM, LCDM, OCDM, TCDM, CHDM) are plotted together with observed data from Kim, et. 
al [9] (open squares, z = 2.85) and Kirkman & Tytler [10] (filled triangles, z = 2.7) (bottom). 
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FIGURE 2. The evolution of the median of the b parameter distribution with redshift for 
SCDM, LCDM, OCDM, & TCDM compared to the observed data by Kim, et. al [9] for lines 
with N HI in the ranges 10 138 - 10 16 cm" 2 (top) and 10 13 1 - 10 14 cm" 2 (bottom). 
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FIGURE 3. Optical depth probability distribution at z = 3. Results for SCDM with varying 
power (agh-i) (top); results for SCDM, LCDM, OCDM, TCDM, & CHDM (bottom). 



spatial resolution. Fig. 3 shows that the distribution narrows for models with lower 
power and that the shape of the distribution varies significantly for the different 
models over the range 0.05 < r < 4 accessible to observations. Comparison of high 
quality observations with high resolution simulations using an ensemble of such 
statistics may soon clarify the physical properties of the intergalactic medium at 
intermediate redshifts when galaxies were young. 
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